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We investigate gravitino dark matter scenarios in which the primordial ^Li production is catalyzed 
by bound-state formation of long-lived negatively charged particles X~ with *He. In the constrained 
minimal supersymmetric Standard Model (CMSSM) with the stau Tj~ as the X~ , the observationally 
inferred bound on the primordial ®Li abundance allows us to derive a rigid lower limit on the gaugino 
mass parameter for a standard cosmological history. This limit can have severe implications for 
supersymmetry searches at the Large Hadron Collider and for the reheating temperature after 
inflation. 



INTRODUCTION 

Big Bang Nucleosynthesis (BBN) is a powerful tool to 
test physics beyond the Standard Model. Recently, it has 
been realized that the presence of heavy long-lived nega- 
tively charged particles X~ can have a substantial impact 
on the primordial light element abundances via bound- 
state formation [1-9]. In particular, when X~ and "^He 
form Coulomb bound states, (^HeX~), too much ^Li can 
be produced via the catalyzed BBN (CBBN) reaction [1] 

i^ReX-) +B ^^^U + X-. (1) 

The formation of {^HeX^) and hence the CBBN produc- 
tion of ^Li becomes efficient at temperatures T ^ 10 keV, 
i.e., at cosmic times t > 10^ s at which standard BBN 
(SBBN) processes are already frozen out. The observa- 
tionally inferred bound on the primordial ^Li abundance 
then restricts severely the X^ abundance at such times. 

A long-lived X~ may be realized if the gravitino is 
the lightest supersymmetric particle (LSP). In particu- 
lar, it is reasonable to consider gravitino LSP scenarios 
within the constrained minimal supersymmetric Stan- 
dard Model (CMSSM) [4, 10-13] in which the gaugino 
masses, the scalar masses, and the trilinear scalar cou- 
plings are parameterized by their respective universal 
values mi/2, mo, and Aq at the scale of grand unifica- 
tion Mgut =i 2 X lO^*^ GeV. Within this framework, the 
lighter stau ti is the lightest Standard Model superpart- 
ncr in a large region of the parameter space and thus a 
well-motivated candidate for the ncxt-to-lightest super- 
symmetric particle (NLSP). Since its couplings to the 
gravitino LSP are suppressed by the (reduced) Planck 
scale, Mp = 2.4 x 10^^ GeV, the stau wiU typically be 
long-lived for conserved R-parity^ and thus -rf can play 
the role of X~ . 

In scenarios with conserved R-parity, the gravitino LSP 
is stable and a promisig dark matter candidate. Graviti- 
nos can be produced efhciently in thermal scattering of 
particles in the primordial plasma. If the Universe, after 



inflation, enters the radiation dominated epoch with a 
high reheating temperature Tr, the resulting gravitino 
density il-^ will contribute substantially to the dark 
matter density ildm [15-17]. 

In this work we calculate the amount of ^Li produced 
in (1) by following the treatment of Ref. [18]. In parti- 
cluar, we employ a recent state-of-the-art result for the 
CBBN reaction cross reaction [5]. The obtained upper 
limit on the X~ abundance from possible ^Li overpro- 
duction vanishes for sufficiently short tx- ■ This allows 
us to extract a lower limit on the universal gaugino mass 
parameter TOi/2 within minimal supergravity scenarios 
where the gravitino is the LSP and the X^ is the Tj~ 
NLSP.^ This limit leads directly to an upper bound on 
Tr since fl^ cannot exceed the observed dark matter 
density. The bounds on mi/2 and Tr derived below de- 
pend on the gravitino mass but are independent of the 
CMSSM parameters. 

Before proceeding, let us comment on the present sta- 
tus of BBN constraints on gravitino dark matter scenar- 
ios with a long-lived charged slepton NLSP. In a recent 
ambitious study [9] it is argued that bound-state forma- 
tion of X^ with protons at T ^ 1 keV might well re- 
process large fractions of the previously synthesized ^Li. 
This seems to relax the bound on the X~ abundance for 
Tx- > 10^ s. However, at present, the uncertainties in 
the relevant nuclear reaction rates in [9] make it difficult 
to decide whether a new cosmologically allowed region 
will open up. In this work we assume that this is not the 
case, in particular, since the '^He/D constraint on elec- 
tromagnetic energy release [19] becomes severe in this re- 
gion and excludes stau lifetimes r^^ > 10^ s [4, 7, 9, 11]. 
Then only the constraint from hadronic energy release on 
D [11, 20-23] can be slightly more severe than the one 
from catalyzed ^Li production [4, 7, 13, 24]. We neglect 
the D constraint in this work since it can only tighten the 
bounds on mi/2 ^-nd Tr, as can be seen, e.g., in Figs. 4 (b- 
d) and 5 of Ref. [13]. For deriving conservative bounds on 



For the case of broken R-parity, see, e.g. [14] 



^ In this work we assume a standard cosmological history with a 
reheating temperature Tr that exceeds the freeze-out tempera- 
ture Tf of the Ti NLSP. 
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mi/2 and Tr , it is thus sufficient to consider the CBBN 
reaction (1) exclusively. 

CATALYZED ^LI PRODUCTION 

The following set of Boltzmann equations [18] describe 
the time evolution of bound-state (BS) formation of X~ 
with ^He and the associated evolution of the primordial 
light elements involved: 



dt 
dYx- 
dt 

dt 

dt 

dYu 
dt 



= {u,v)sYs - Tx-Yjis - {cjcv)sY^sYt,, (2a) 

= -{a,v)sYs - Tx-Yx- + {acv)sYj,sYT,, (2b) 

= -{a,v)sYs+Vx-Y^s, (2c) 

= {acv)sY^sYT,, (2d) 

= -{gcv)sY^sYt,. (2e) 



Here wc scale out the expansion of the Universe by defin- 
ing the yield Y^ = iii/s where rzj is the number density 
of species i and s = 27r^ 5f*s is the entropy density. 

In particular, IbS; Yx- , Yi^^, leu, and Id denote the 
yields of the (^HeX~) bound state, free , free ^He, ^Li 
produced in CBBN, and D, respectively. The quantity 
Ys = (Yx-Y4ii^ — 188^7) parameterizes the competition 
between recombination and photo^iissociation of bound 
states. For the latter, one defines Y^ = rij/s with [2] 



Vile 



3/2 



(3) 



2C(3) V27rTy' ^ 

where E^, = 337.33 kcV [5] is the (^HeX") binding energy 
and = 2^(3)T'^/7r^. Furthermore, Tx- = t^- denotes 
the total decay width of X~ . 

The CBBN reaction cross section for the process (1) 
has recently been computed with an advanced method 
from nuclear physics [5] 



(acv) = 2.37 X 10** (1 - O.34T9) Tg^'e 



2/3^-5.33T- 



(4) 



which is given in imits of iV^"'^cm^s~^mole~'^ with Tg de- 
noting the temperature in units of 10^ K. The recombi- 
nation cross section of X~ with ^He is estimated as [2] 



(cTrW) 



3e^ m\\fm^ 
with ma = 3.73 GeV [25] and Z^, = 2? 



(5) 



^ Equation (5) assumes a radiative capture of X~ into the IS 
bound state of a point-like ot particle. We use, however, 
obtained numerically in [5] rather than the Bohr-like formula 
Eg ~ Z2Q2ma/2 = 397 keV; ma < ra^- ■ 




FIG. 1: Contour-lines of ^Li/H produced in CBBN obtained 

by solving (2) (solid lines) and by using the Saha type approx- 
imation for Ybs instead of computing (2a) (dashed lines). 



We solve (2) using as initial conditions the respec- 
tive X~ yield prior to decay, Y^^ , and the SBBN 
output values of the computer code PArthENoPE [26]: 
yp = 4n4Hc/nb = 0.248, D/H = 2.6 x 10-^ ^Li/H = 
1.14 X 10"^**, and n-pln\j = 0.75; furthermore, 5* = 3.36 
and g^:S = 3.91. While the variation of the D and *He 
abundances from their SBBN values are negligible, the 
catalyzed fusion of ^Li is substantial as shown in Fig. 1 
by the contour-lines of ^Li/H = Yej^^s/np (solid lines). 
Contrasting with the observationally inferred upper limit 
on the primordial ^Li abundance [27], 



=^Li/H|obs<2xlO-", 



(6) 



one sees clearly that ^Li/H]cBBN can be far in excess. 

The dashed lines in Fig. 1 show the solution of (2) 

where instead of (2a) the Saha type equation Ybs = 
Yi-HeYx- /Yj is used as an approximation for the bound- 
state abundance. The obtained overestimation of the 
*'Li abimdance demonstrates the importance of the use 
of the Boltzmann equation (2a). However, focusing on 
we will read off the relevant constraint 
in the region tx- < 10^ s in which the slope of the ^Li 
contours is very steep. Therefore, the use of (2a) instead 
of the Saha type equation is an improvement on the con- 
ceptional side which has only a marginal effect on the 
bounds to be derived. By the same token, and from the 
comparison of our results with [4], we find that also the 
destruction of ^Li due to X~ decays affects those bounds 
only marginally. 
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LOWER LIMIT ON M1/2 

Applying the above results to gravitino dark mat- 
ter scenarios with the lighter stau ti as the NLSP, we 
now derive the conservative lower limit on mi/2- The 
stau NLSP with a mass of m^i decouples from the pri- 
mordial plasma with a typical yield of Y-^'^ ^ 7 x 
10-"(m^=J100 GeV) [28]. With = Y^^" 12, we find 
from Fig. 1 that the amount of ^Li produced in CBBN 
can be in agreement with (6) only for stau lifetimes of 



= Tx- ^ 5 X 10^ s. 
As can be seen from the supergravity prediction 



(7) 



m| 
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[GeV] 



1000 



the requirement (7) implies a lower limit on the splitting 
between my^ and provided niTf^ < 0(1 TeV). Be- 
cause of this hierarchy, the factor (1 — Ti^/rTi-^)^"* can 
be neglected in the following. 

Let us now turn to the CMSSM. In the region in 
which Ti is the NLSP, we find 



(9) 



by scanning over the following parameter range: 

TOi/2 = 0.1-6 TeV, 
tan /? = 2-60, 
sgn^ = ±1, 
— 4too < Aq < 4too 

with mo as large as viable for a t\ NLSP. Here tan /3 is 
the ratio of the two MSSM Higgs doublet vacuum expec- 
tation values and /i the Higgsino mass parameter.'* 

For small left-right mixing, t\ ~ tr , (9) can be under- 
stood qualitatively from the estimate for the mass of the 
right-handed stau TOtr near the electroweak scale [30] 



^.\hra\ 



/2 



sin Qy/f^i cos 2/3 



(10) 



since rn^ ^ rr? 



Q ^ "'1/2 ^ large part of the ri NLSP region. 
In fact, (9) tends to be saturated for larger mo, i.e., in 
the stau-neutralino-coannihilation region where the mass 
of the lightest neutralino m^^o ~ m^i ■ This can be un- 
derstood since the neutralino is bino-like in this region 



FIG. 2: The shaded region indicates cosmologically disfavored 
mi/2 values. Below the dashed line, mg > rrir^ is possible. 



so that 



0.18m2/2.^ 



In the remaining part of the 



stau NLSP region, smaller values of m?^ satisfying, e.g., 



0.15m^y2 '^^'^ easily be found. 



To be on the conservative side, we set the stau NLSP 
mass m^ij to its maximum value at which (9) is saturated: 
m~ = 0.21mJ^2- Then, constraint (7) together with (8) 
yields 



"^•1/2 



> 0.9 TeV 



10 Gev) 



2/5 



(11) 



which is shown in Fig. 2. The shaded region is disfavored 
by (6). Below the dashed line, m^ > m^r^ is possible. 

Since for a ri NLSP typically m§ <C '7i^,2j it is the 
gaugino mass parameter mi/2 which sets the scale for 
the low energy superparticle spectrum. Thus, depending 
on m^, the bound (11) implies rather high values of the 
superparticle masses. This is particularly true for the 
masses of the squarks and the gluino since their renor- 
malization group running from Mqut to Q ~ 0(1 TeV) 
is dominated by M3((5) — mi/2as((9)/as(-^^GUT)- Since 
these masses govern the size of the total cross section 
for the production of superparticles at the Large Hadron 
Collider (LHC), the cosmologically favored region for 
"^G ^ 1*-* GeV is associated with a mass range that will 
be very diSicult to probe at the LHC. 

Let us stress at this point that the bounds (7) and (11) 
and their severe implications for phenomenology at the 



* We employ SPheno 2.2.3 [29] to compute the low energy mass 
spectrum using rrtt = 172.5 GeV for the top quark mass. In 
addition, we use the Standard Model parameters rrtij(mb)'^^ = 
4.2 GeV, a?^(mz) = 0.1172, (m^) = 127.932. 



This estimate is relatively independent of tan /3 and valid in the 
mi/2 region in which also the LEP bound on the Higgs mass [25], 
rrih > 114.4 GeV, is respected. 
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LHC arc valid only for the assumed standard cosmolog- 
ical history with Tr > Tf and the associated consid- 
ered values of Y-^'^'^ ^ 7 x 10~^^. For example, non- 
standard entropy production after the thermal ti NLSP 
freeze out and before BBN might dilute the stau abun- 
dance prior to decay. Thereby, the mi/2 limit can be 
relaxed [13]. Also for the case of inflation with a low re- 
heating temperature, Tr < Tf, one can obtain a stau 
abundance prior to decay that respects the ^Li con- 
straint even for t^-^ ^ 5 x 10^ s [18]. Thus, for a 
non-standard cosmolocial history, an observation of staus 
with T^Tj > 5 X 10^ s and other CMSSM phenomenology 
at the LHC could still be viable even in gravitino LSP 
scenarios with > 10 GeV. 



Note that the constraint (14) relies on thermal grav- 
itino production only. In addition, gravitinos are pro- 
duced in stau NLSP decays with the respective density 

nNTP^2 ^ rn^Y^^;'s{To)hyp, , (15) 

where pc/[s{To)h^] = 3.6 x IQ-^ GeV [25]. While the 
precise value of Y-^'^ depends on the concrete choice of 
the CMSSM parameters, the upper limit (14) can only 
become more stringent by taking into account. For 

exemplary CMSSM scenarios, this can be seen from the 
{mi/2,nio) planes shown in Figs. 4 and 5 of Ref. [13]. 
These figures illustrate that the severe limits (11) and 
(14) are conservative bounds. 



UPPER BOUND ON Tr 

The amoimt of gravitinos produced in thermal scatter- 
ing is sensitive to the reheating temperature Tr and to 
the masses of the gauginos and hence to mi/2 [16]. For 
a standard cosmological history, the associated gravitino 
density can be approximated by^ 

This follows from Eq. (3) of Ref. [16]. Here we use that 
the running gaugino masses Mj associated with the gauge 

groups SU(3)c, SU(2)l, and U(1)y satisfy M3 : M2 : 
Ml ~ 3 : 1.6 : 1 at a representative scale of 10^ GeV at 
which we also evaluate the respective gauge couplings. 
Furthermore, we only need to take into account the pro- 
duction of the spin- 1/2 components of the gravitino since 
(11) implies M2/3m| > 1 for > 1 GeV. 

For a given the reheating temperature Tr is lim- 

ited from above because cannot exceed th(^ dark 

Cr 

matter density [25] ^I'^h'^ = ClOSlg where h is the 
Hubble constant in units of 100 km Mpc~^s~^ . Requiring 

i^o'h'^ < 0.126 (13) 

and using the derived lower bound (11) allows us to ex- 
tract the conservative upper limit: 

r„<4.9xlO'Gev(j^)"\ (14) 

This constraint is a slowly varying function of mg: 
{m^/lOGeVy/^ = 0.6 - 2.5 for = 1 GeV - 1 TeV. 
Therefore, (14) poses a strong bound on Tr for the natu- 
ral gravitino LSP mass range in gravity-mediated super- 
symmetry breaking scenarios. 



CONCLUSION 

We have considered the catalysis of ^Li production 
in CMSSM scenarios with the gravitino LSP and the 
stau NLSP. Within a standard cosmological history, the 
calculated ^Li abundance drops below the observational 
limit on primordial ^Li for t^^ ^ 5 x 10^ s. Taken at 
face value, we find that this constraint translates into 
a lower limit TOi/2 > 0.9 TeV(TOg/10 GeV)^/^ in the 
entire natural region of the CMSSM parameter space. 
This implies a conservative upper bound Tr < 4.9 x 
10^GeV(mg/10GeV)i/5. The bounds on toi/2 and Tr 
not only confirm our previous findings [13] but are also 
independent of the particular values of the CMSSM pa- 
rameters for the considered n NLSP abundances. 

We are grateful to T. Plehn, S. Reinartz, and A. Weber 
for valuable discussions. 

Note added: After submission of this work, a substan- 
tially revised version (v3) of [9] together with [33] ap- 
peared on the arXiv. The results of these works affect our 
limits only mildly. Because of the huge effect of (1) on the 
^Li abundance, our relatively simple treatment of CBBN 
is sufficient for our purposes. This is also confirmed by a 
direct comparison of our data with Figs. 1 and 2 of the 
more elaborate CBBN treatment in [33] for Sh < 3 x 10"^ 
{rur^ < 2.7 TeV, i.e., mi/2 < 6 TeV) [23] at the relevant 
times of t ~ few x 10"^ s. For a given '^Li/Hjobg bound, 
the effect on the r^^ limit is less than a factor of 1.5. In 
addition, adopting ^Li/HJobs < 4x 10"" (2.7 x 10"!°) as 
used in [33], the numbers in our Eqs. (7), (11), and (14) 
change respectively to 6 x 10^ (10^), 0.87 (0.78), and 
5.3 X 10^ (6.5 X lO''). Furthermore, by taking into ac- 
count the uncertainties in the relevant nuclear reaction 
rates, it is shown explicitly in Fig. 14 in v3 of [9] and 
in Fig. 5 in [33] that cosmologically allowed regions for 
> 10^ s are indeed extremely unlikely (< 1%) for 



For a discussion on the definition of Tj^, see Sec. 2 in Ref. [13]. 



Gravitino production from inflaton decay can also be substantial; 
see, e.g., [31, 32]. This can further tighten the bound (14). 
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ydec > 7 X 10-"(m?^/100 GeV) even with /o^ as small 
as 3 X 10~^. Only with a finely tuned mf^-niQ degener- 
acy leading to Sh — > and /em — > can any bound on 
energy release and, in particular, the one from ^He/D be 
evaded. 
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